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Abstract 

Background: The genus Haematopinus contains 21 species of blood-sucl<ing lice, parasitizing both even-toed ungulates 
(pigs, cattle, buffalo, antelopes, camels and deer) and odd-toed ungulates (horses, donkeys and zebras). The mitochondrial 
genomes of the domestic pig louse, Haematopinus suis, and the wild pig louse, Haematopinus apri, have been sequenced 
recently; both lice have fragmented mitochondrial genomes with 37 genes on nine minichromosomes. To understand 
whether the composition of mitochondrial minichromosomes and the gene content and gene arrangement of 
each minichromosome are stable within the genus, we sequenced the mitochondrial genome of the horse 
louse, Haematopinus asini. 

Methods: We used a PCR-based strategy to amplify four mitochondrial minichromosomes in near full-length, 
and then amplify the entire coding regions of all of the nine mitochondrial minichromosomes of the horse 
louse. These amplicons were sequenced with an lllumina Hiseq platform. 

Results: We identified all of the 37 mitochondrial genes typical of bilateral animals in the horse louse, 
Haematopinus asini; these genes are on nine circular minichromosomes. Each minichromosome is 3.5-5.0 kb in 
size and consists of a coding region and a non-coding region except R-nad4L-rrnS-C minichromosome, which 
contains two coding regions and two non-coding regions. Six of the nine minichromosomes of the horse louse 
have their counterparts in the pig lice with the same gene content and gene arrangement. However, the gene 
content and arrangement of the other three minichromosomes of the horse louse, including R-nad4L-rrnS-C, 
are different from that of the other three minichromosomes of the pig lice. 

Conclusions: Comparison between the horse louse and the pig lice revealed variation in the composition 
of mitochondrial minichromosomes within the genus Haematopinus, which can be accounted for by gene 
translocation events between minichromosomes. The current study indicates that inter-minichromosome 
recombination plays a major role in generating the variation in the composition of mitochondrial minichromosomes 
and provides novel insights into the evolution of fragmented mitochondrial genomes in the blood-sucking lice. 
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Background 

Haematopinus (Leach, 1815) is the only genus in the family 
Haematopinidae of the suborder Anoplura, known as the 
blood-sucking lice [1-3]. There are 21 species in the genus 
Haematopinus, of which 19 species parasitize even- toed 
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ungulates (order Artiodactyla) such as pigs, cattle, buffalo, 
antelopes, camels and deer, whereas the other two species 
parasitize odd-toed ungulates (order Perissodactyla) 
such as horses, donkeys and zebras [1,4]. Haematopi- 
nus species are vectors of several severe infectious dis- 
eases in rural tropical areas such as African swine fever 
[5-7], swinepox [8], hog cholera and eperythrozoonosis 
[9,10], and anaplasmosis [11]. 
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The typical mitochondrial (mt) genome organization of 
insects and other bilateral animals consists of a single circu- 
lar chromosome, 13-20 kb in size, with 37 genes and a 
control region [12,13]. Fragmented mt genomes, however, 
have been found in five species of blood-sucking lice: hu- 
man body louse, Pediculus humanus, human head louse, 
Pe. capitis, human pubic louse, Pthirus pubis, domestic pig 
louse, Haematopinus suis and the wild pig louse, H, apri 
[14-16]. The extent of mt genome fragmentation, the num- 
ber of minichromosomes and the distribution of genes on 
the minichromosomes vary remarkably between genera of 
the blood-sucking lice, but are the same for the species 
within the same genus [14-16]. For instance, both the hu- 
man body louse, Pe, humanus, and the human head louse, 
Pe, capitis, have 20 mt minichromosomes and an identical 
pattern for the distribution of mt genes on these minichro- 
mosomes [14,16]. The domestic pig louse, H, suis, and the 
wild pig louse, H, apri, have nine mt minichromosomes 
and an identical pattern for the distribution of mt genes too 
[15]. The human head louse and the human body louse, 
however, are very closely related and had their most recent 
common ancestor (MRCA) -107,000 years ago [17,18]; so 
are the domestic pig louse and the wild pig louse, which 
had their MRCA --9,000 years ago [19]. 

To understand whether the composition of mt mini- 
chromosomes, and the gene content and gene arrange- 
ment in each minichromosome are indeed conserved 
among species of the same genus, we sequenced the mt 
genome of the horse louse, Haematopinus asini, and 
compared it with those of the pig lice, H, suis and H, 
apri, H, asini parasitizes horses {Equus caballus) and 
two other odd-toed ungulates: donkeys (£. asinus) and 
plains zebras (£. burchelli) [20]. We found that the horse 
louse differs from the pig lice in the distribution of mt 
genes on three of the nine minichromosomes and the 
intra-genus variation can be explained by gene transloca- 
tion between minichromosomes. 

Methods 

Sample collection, DNA extraction, PGR amplification and 
DNA sequencing 

Specimens of H, asini were collected from horses at Kemps 
Creek, Sydney, AustraUa in 2008 (sample # B2448). 
The study does not involve the ethical treatment of 
horse, other than collecting lice from their body sur- 
face. Genomic DNA was extracted from individual 
louse specimens with DNeasy Tissue kit (QIAGEN). Four 
pairs of primers, 12SA-12SB, 16SF-Lxl6SR, mtd6-mtdll 
and mtdl6-mtdl8, were used to amplify fragments of rrnS 
(350 bp), rrnL (300 bp), coxl (600 bp) and cox2 (230 bp) 
(see Additional file 1). These fragments were sequenced 
with AB3730xl sequencers. Four pairs of outbound primers 
(forward and reverse), 12sB2448F-12sB2448R, 16sB2448F- 
16sB2448R, coxlB2448F-coxlB2448R and cox2B2448F- 



cox2B2448R (see Additional file 1), were designed from the 
sequences of the rrnS, rrnL, coxl and cox2 fragments, re- 
spectively. PCRs with these specific primers amplified four 
mt minichromosomes in near fijill-length that contained 
each of the four genes, 3.5 kb, 3.8 kb, 5.0 kb and 4.1 kb in 
size respectively (Figure lA). These amplicons were 
sequenced partially with AB3730xl sequencers. Se- 
quences from the non-coding regions of these four 
minichromosomes were aligned with Clustal X [21]. A 
forward primer B2448F and a reverse primer B2448R 
were designed from conserved motifs in the non-coding 
regions that flank the coding regions of the four minichro- 
mosomes above (see Additional file 1). A mixture of PGR 
ampUcons ranging from 350 bp to 2,900 bp was ob- 
tained with the primer pair B2448F-B2448R; these 
amplicons were expected from the coding regions of all mt 
minichromosomes of the horse louse (Figure IB). Ampli- 
cons generated with B2448F-B2448R and those with pri- 
mer pairs 12sB2448F-12sB2448R, 16sB2448F-16sB2448R, 
coxlB2448F-coxlB2448R and cox2B2448F-cox2B2448R 
were sequenced with Illumina Hiseq 2000 platform at the 
BGI, Hong Kong. 

Takara Ex Taq was used in the initial short PGRs with 
the following cycling conditions: 94°G for 1 min; 35 cy- 
cles of 98°G for 10 sec, 45°G for 30 sec, 72°G for 1 min; 
and a final extension of 72°G for 2 min. TaKaRa LA Taq 
was used in the long PGRs with the cycling conditions: 
94°G for 1 min, 35 cycles of 98°G for 10 sec, 55-65°G 
(depending on primers) for 40 sec, 68°G for 4 min; and 
72°G for 8 min. Negative controls were executed with 
each PGR experiment. PGR amplicons were checked by 
agarose gel (1%) electrophoresis; the sizes of PGR 
amplicons were estimated by comparison with molecu- 
lar markers. Wizard SV Gel and PGR Glean-up System 
(Promega) were used to purify PGR amplicons for 
sequencing. 

Assembly of Illumina sequence-reads, gene identification 
and verification of mitochondrial minichromosomes 

lUumina sequence-reads were assembled de novo with 
Geneious (Version 6.1.6, Biomatters); the assembly pa- 
rameters were 98% and 50 bp for minimum overlap 
identity and minimum overlap, respectively. tRNA genes 
were identified using program tRNAscan-SE [22] and 
ARWEN [23]. Protein-coding genes and rRNA genes were 
identified with BLAST searches of GenBank [24,25]. Identi- 
cal sequences shared between genes were searched with 
program Wordmatch [26]. The length of identical se- 
quences shared by chance between genes was assessed by 
analyzing randomly-extracted, unrelated DNA sequences 
from GenBank of the same sizes as our experimental se- 
quences. Take tRNA genes, which are ~ 70 nt in length, as 
an example. We extract, randomly, a set of unrelated DNA 
sequences (n > 50) from GenBank; each of these sequences 
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Figure 1 PGR amplicons from the mitochondrial genome of the horse louse, Haematopinus asini. (A) Amplicons generated with tine liorse-louse- 
specific primers, 12sB2448F-12sB2448R (lane 2), 16sB2448F-16sB2448R (lane 3), coxl B2448F-coxl B2448R (lane 4) and cox2B2448F-cox2B2448R (lane 5) 
from four mitochondrial minichromosomes. Lane 1 and lane 6: 100-bp Ladder and 1-kb Ladder (Biosciences). (B) Amplicons generated with the primer 
pair B2448F-B2448R from the coding regions of all of the mitochondrial minichromosomes of the horse louse (lane 2). Lane 1: 500-bp DNA Ladder 
(Tiangen). (C) PCR verification of the mt minichromosomes of the horse louse. Lane 1 and 12: 100-bp ladder. Lane 2 and 13: 1-kb ladder. 
Lane 3-1 1: PCR amplicons from the nine minichromosomes of the horse louse: K-nad4-atp8-otp6-N, nad2-l-coxl-L2, D-Y-cox2-S]-S2-P-cox3-A, 
E-cob-V, Q-nadl -T-G-nad3-W, H-nad5-F-nod6, M, LrrrnL and R-nad4L-rrnS-C. Genes from which PCR primers were designed are in bold. 



is 70 nt. We then run these sequences in Wordmatch and 
identify the size of the longest identical sequences shared 
between any pairs of the extracted sequences. We work out 
the average (~7 nt in the case of tRNA genes) and use it as 
an indication of chance expectation of the length of identi- 
cal sequences for tRNA genes. 

The size and circular organization of each mt mini- 
chromosome identified by sequence-read assembly were 
verified by PCR. A pair of outbound primers (forward 
and reverse) was designed from the coding region of 
each minichromosome (see Additional file 2). PCRs with 
these primers amplified each minichromosome in full- 
length or near full-length if the minichromosomes had a 
circular organization. The amplicons generated with pri- 
mer pairs nad4f-nad4r, cobf-cobr, nad3f-nad3r, nadSf- 
nadSr and metf-metr were also sequenced with Illumina 
Hiseq 2000 platform. PCR set-up, cycling conditions, 
agarose gel electrophoresis and size measurement were 
the same as the long PCRs described above. 

Results and Discussion 

Mitochondrial genome of the horse louse, Haematopinus 
asini 

We obtained 897,097 Illumina sequence-reads (pair end, 
180-bp inserts) from the amplicons of the mt genome of 
the horse louse, H. asini (Table 1). The sequence-reads are 
all 90 bp each in length. We assembled these sequence- 
reads into contigs and identified all of the 37 mt genes 
typical of bilateral animals in H, asini, distributed on nine 
circular minichromosomes (Figure 2; Figure IC). The mt 
minichromosomes of the horse louse are 3.5-5.0 kb in size 
(Figure IC) and are the largest among those of the sucking 
lice known, due to their expanded non-coding regions (see 
below). Each minichromosome of the horse louse consists 
of a coding region and a non-coding region except the 



R-nad4L-rrnS-C minichromosome, which has two coding 
regions and two non-coding regions (Figure 2). There are 
1-8 genes in each coding region, varying in size from 66 bp 
for trnM minichromosome to 2,699 bp for nad2-trnI-coxl- 
trnL2 minichromosome (Table 1). With the exception of 
trnT, nadl and trnQ, all of the mt genes have the same 
orientation of transcription relative to the non-coding 
region (Figure 2). The nucleotide sequences of the mt 
minichromosomes of H, asini were deposited in GenBank 
under accession numbers KF939318, KF939322, KF939324, 
KF939326 and KJ434034-KJ434038 (Table 1). 

We sequenced the non-coding regions of all of the 
nine mt minichromosomes of the horse louse in full 
length, which range from 2,005 bp to 3,264 bp (Figure 3). 
The horse louse is the first species of sucking lice for 
which the full-length non-coding regions of all mt 

Table 1 Mitochondrial minichromosomes of the horse 
louse, Haematopinus asini, identified by Illumina 
sequencing 

Minichromosome GenBank Size of Number of Mean 





accession 
number 


coding 
region (bp) 


Illumina 
sequence- 
reads 


coverage 


D-y-cox2-SrS2-P-cox3-A 


KF939318 


1871 


150403 


3420 


E-cob-V 


KJ434038 


1218 


281565 


6435 


H-nad5-F-nad6 


KJ434037 


2229 


84549 


1579 


K-nad4-atp8-atp6-N 


KJ434035 


2282 


1009104 


18660 


LrrrnL 


KF939322 


1289 


120832 


7590 


M 


KJ434036 


66 


409434 


10790 


nod2-l-coxl-L2 


KF939324 


2699 


133014 


2990 


Q-nodl-T-G-nod3-W 


KJ434034 


1505 


487487 


8846 


R-nad4L-rmS-C 


KF939326 


1142 


16,059 


385 


Total 




14301 


2692447 


60695 
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Figure 2 Mitochondrial genome of the horse louse, Haematopinus asini. The name, transcription orientation and lengtli (bp) of eacli gene 
are indicated. Non-coding regions are in blacl<. Abbreviations of gene names are: coxl-3 for cytoclirome c oxidase subunits 1-3; cob for cytoclirome b; 
nadl-S and nad4L for NADH deliydrogenase subunits 1 -5 and 4 L; and rrnS and rrnL for small and large ribosome RNA subunits. tRNA genes 
are labeled with the single-letter abbreviations of their corresponding amino acids. Numbers indicate the length of each corresponding 
gene. Minichromosomes shown with asterisk symbols (*) have different gene content and gene arrangement compared with the pig lice, 
Haematopinus suis and Haematopinus apri [15]. 



minichromosomes were sequenced. The horse louse 
has the longest non-coding regions among the sucking lice 
known; previously the longest non-coding region was 
2,370 bp, noted in the pig lice [15]. As in the human lice 
and the pig lice, there is an AT-rich motif (45 bp, 100% A 
and T) in the non-coding region upstream the 5 '-end of 
coding region and a GC-rich motif (78 bp, 60% C and G) 
downstream the 3 '-end of the coding region (Figure 3). 
The size variation among the nine non-coding regions of 
the horse louse is due to size variation in the section 
upstream the coding region from the AT-rich motif to the 
primer B2448F (Figure 3). Excluding this section, the non- 
coding regions of the minichromosomes have ~ 96% pair- 
wise identity to each other. 

Shared identical sequences and recombination between 
mt genes in the horse louse 

Ten pairs of mitochondrial genes share stretches of 
identical sequences longer than expected by chance in 



the blood-sucking lice of humans and pigs, providing 
unequivocal evidence for DNA recombination between 
mt genes and between minichromosomes in these lice 
[14,16]. We found that nine pairs of mt genes in the 
horse louse, H. asini, also share stretches of identical se- 
quences longer than expected by chance (Table 2). trnLj 
and trnL2 are the only pair of genes that share longer 
than expected identical sequences in all of the three spe- 
cies of human lice and the two species of pig lice, indi- 
cating recombination between these two genes has been 
maintained since their MRCA, ~ 65 million years ago 
(Mya) [27]. This is also the case for the horse louse al- 
though the 15-bp identical sequence shared between 
these two genes is the shortest among the sucking lice 
known. Any two non-homologous tRNA genes might be 
expected to share identical sequences, ~ 7 bp in size, by 
chance; the identical sequence shared between trnL j and 
trnL2 in the horse louse is twice as long as expected by 
chance (Table 2; see Additional file 3). As in the two pig 
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Figure 3 Alignment of the full-length non-coding regions of nine mitochondrial minichromosomes of the horse louse, Haematopinus 
asini. B2448F and B2448R are the primers used to amplify tine entire coding regions of all mitochondrial minichromosomes of the horse louse. 
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Table 2 The longest stretches of identical sequence shared between mitochondrial genes in the horse louse, 
Haematopinus asini 

The longest stretches of identical sequence shared (bp) 



Pairs of genes 


Horse louse 


Pig lice 




Human lice 




Animals with typical mt genome organization 


Haas 


Has 


Haap 


Peh 


Pec 


Ptp 


Bm 


Cb 


Hm 


Dy 


cob 


coxl 


21 


11 


11 


10 


10 


10 


12 


13 


14 


13 


nodi 


nad4 


20 


9 


12 


14 


14 


N/A 


12 


16 


13 


12 


nod4 


cob 


20 


10 


10 


11 


11 


N/A 


10 


12 


12 


12 


nod2 


nods 


25 


10 


10 


14 


14 


10 


12 


16 


12 


15 


nod2 


nods 


19 


12 


12 


10 


10 


10 


10 


12 


11 


11 


rrnL 


rrnS 


10, 33 


10 


9 


11 


11 


10 


33 


10 


15 


16 


rrnL 


trnN 


15 


9 


9 


9 


9 


N/A 


10 


8 


9 


8 


trnL^ 


trnL2 


15 


16,10,9 


16,10,9 


33,32 


33,32 


35,32 


7 


6 


7 


10 


tmJ 


trnP 


27 


26 


26 


7 


7 


N/A 


6 


8 


8 


9 



Note: Abbreviations of species names are: Haos, Haematopinus asini (horse louse); Has, Haematopinus suis (domestic pig louse); Haap, Haematopinus apri (wild pig 
louse); Peh, Pediculus humanus (human body louse); Pec, Pediculus capitis (human head louse); Ptp, Pthirus pubis (human pubic louse); Bm, Bothriometopus 
macrocnemis (screamer louse); Cb, Campanulotes bidentatus (pigeon louse); Hm, Heterodoxus macropus (wallaby louse); Dy, Drosophila yakuba (fruitfly); N/A, Not 
available, mt. Mitochondrial; Stretches of shared identical sequences longer than expected by chance are in bold. 



lice, trnTdind trnP in the horse louse share 27-bp identi- 
cal sequence, which is four times as long as expected by 
chance (Table 2; see Additional file 3). In the three hu- 
man lice and other animals, trnT zxvd trnP share identi- 
cal sequences, 6-9 bp long, which are expected by 
chance. Recombination between trnT and trnP, thus, is 
likely a derived feature for the genus Haematopinus, 
Seven other pairs of mt genes also share identical se- 
quences 1.5-3 times longer than expected by chance in 
the horse louse but not in the pig lice nor in the human 
lice (Table 2). Recombination between these seven pairs 
of genes, therefore, is likely a derived feature for the 
horse louse only (Table 2). 



Variation in mt minichromosome composition between 
the horse louse and the pig lice 

Six of the nine minichromosomes of the horse louse, H. 
asini, have the same gene content and gene arrangement 
as their counterparts of the pig lice, H, suis and H, apri 
(Figure 2) [15]. These minichromosomes are apparently 
ancestral to Haematopinus species and thus have been 
retained in both the horse louse and the pig lice. The 
other three minichromosomes of the horse louse, how- 
ever, are not present in the pig lice (Figure 2). In the pig 
lice, one of the minichromosomes has three genes, trnH- 
nadS-trnF [15]. the horse louse, however, the minichro- 
mosome that has these three genes also has nad6 gene 




Minichromc^omes 
seen in pig lice 



Minrchronnosonries 
seen in horse 
louse 



Figure 4 An inter-minichromosome recombination model that accounts for the variation in the composition of mitochondrial 
minichromosomes between the horse louse, Haematopinus asini, and the pig lice, H. suis and H. apri. Minichromosomes shown 
with asterisk symbols {*) were seen in a Polyplax rat louse [31]. 
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downstream trnF with a gap of 3 bp in between (Figure 2). 
Similarly, another minichromosome of the pig lice has two 
genes, rrnS-trnC, In the horse louse, however, the mini- 
chromosome that has these two genes has trnR-nad4L up- 
stream rrnS with a 417-bp non-coding region in between 
(Figure 2). Furthermore, in the pig lice, trnM gene is on a 
minichromosome with nad6 and trnR-nad4L [15]. In the 
horse louse, however, trnM is alone on its own minichro- 
mosome (Figure 2). 

How was the intra-genus variation in mt minichromosome 
composition generated? 

Our comparison of the mt genomes of the horse louse 
and the pig lice revealed variation in the composition of 
mt minichromosomes within the genus Haematopinus, 
Several previous studies also compared mt genomes be- 
tween species of sucking lice in the same genus. The hu- 
man head louse and the human body louse in the genus 
Pediculus have identical mt minichromosome compos- 
ition [17,18], so are the domestic pig louse and the wild 
pig louse in the genus Haematopinus [19]. The current 
study compared Haematopinus species that infest mam- 
mals distinct from one another. Haematopinus asini par- 
asitizes exclusively horses, donkeys and zebras, which 
are odd- toed ungulates (order Perissodactyla). Haemato- 
pinus suis and H apri, however, parasitize exclusively 
domestic pigs and wild pigs, which are even-toed ungu- 
lates (order Artiodactyla). These two lineages of ungulate 
mammals had their MRCA 63-83 Mya [28-30]. 

A very recent study by Dong et al. showed that two 
species of rat lice in the genus Polyplax also differ in 
the composition of mt minichromosomes [31]. Together, 
these studies indicate that intra-genus variation in mt 
minichromosome composition is likely common in blood- 
sucking lice. Furthermore, these studies provided oppor- 
tunities to look into how fragmented mt genomes evolved 
in the blood-sucking lice. The typical single-chromosome 
mt genomes are highly conserved in genome organization, 
gene content and gene arrangement in the vast majority 
of insects [12,32]. The fragmented mt genomes of the 
sucking lice known to date, however, showed very limited 
conservation in terms of the number of minichromo- 
somes, and the gene content and gene arrangement in 
each minichromosome. The current study indicates that 
inter- minichromosome recombination plays a major role 
in the fast evolution of fragmented mt genomes in the 
blood-sucking lice [33]. The variation between the horse 
louse and the pig lice in mt minichromosome composition 
can be accounted for parsimoniously by two events of 
inter-minichromosome recombination. Firstly, a recombin- 
ation event translocated R-nad4L from a minichromosome 
that contained R-nad4L-nad6-M (a minichromosome seen 
in the pig lice, [15]) to a minichromosome that contained 
rrnS-trnC and generated a minichromosome with two 



coding regions and two non-coding regions in the horse 
louse: R-nad4L in one coding region whereas rrnS-trnC in 
another coding region (Figure 4). The minichromosome 
that contained rrnS-trnC was seen in the Haematopinus pig 
lice [15] and a Polyplax rat louse [31], and thus can be 
inferred to be ancestral to Haematopinus species. Secondly, 
another recombination event translocated nad6 from the 
R-nad4L-nad6-M minichromosome to a minichromosome 
that contained trnH-nadS-trnF, and generated a minichro- 
mosome with both nad6 and trnH-nadS-trnF in the horse 
louse (Figure 4). The minichromosome that contained 
trnH-nadS-trnF W2is seen in the Haematopinus pig lice [15] 
and a Polyplax rat louse [31], and can be inferred to be 
ancestral to Haematopinus species. Finally, the loss of R- 
nad4L and nad6 from R-nad4L-nad6-M minichromosome 
led to a minichromosome with only trnM, seen in the horse 
louse (Figure 4). 

Conclusions 

We sequenced the mt genome of the horse louse, H 
asini. We found all of the 37 mt genes typical of insects 
and other bilateral animals; these genes are on nine cir- 
cular minichromosomes. Each minichromosome is 3.5- 
5.0 kb in size and contains 1-8 genes. Three of the nine 
minichromosomes of the horse louse differ from those 
of the pig lice in gene content and gene arrangement, re- 
vealing variation in the composition of mt minichromo- 
somes among species of the genus Haematopinus, We 
propose that inter-minichromosome recombination can 
cause gene translocations and likely plays a major role in 
generating the variation in the composition of mt mini- 
chromosomes observed in the Haematopinus and other 
blood-sucking lice. 
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